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ABSTRACT 


© 1.0. Introduction 


Seismic attributes are quantitative measures of seismic characteristics of interest (Liner, 2004; Oyeyemi, and 
Aizebeokhai, 2015; Adewoye, et al, 2015). Seismic attributes need to be integrated in some interpretations as they 
can be used as tools for predicting reservoir geometry and possibly displaying lateral changes in thickness including 
fluid contacts (Steiner, and Lockhart, 1988). Seismic attributes analysis involves the procedure used to extract 
corresponding subsurface geological information from seismic sections (Allstair, 2011; Avseth, et al, 2005). 
Structural and stratigraphic traps could be very subtle and are therefore difficult to map accurately (Omoja, and 
Obiekezie, 2019). The geometrical attributes are used for structural and stratigraphic interpretations of seismic data. 
Seismic attributes are widely used for detection of hydrocarbon reservoirs (Pramanik, et al, 2002; Van, 2000; Vig et 


al, 2002). 


Improved methodologies are needed to discriminate between true hydrocarbon indicators and non-indicators by 
imaging the detail subsurface structure with a view to delineate new hydrocarbon zones that might be hidden in 
structural traps for the development programmes of the field. Seismic attributes analysis being an integral part of 
3D seismic interpretation is one of these advancements (Chopra, and Marfurt, 2005). Several researchers have made 
enormous contributions based on seismic attributes for enhanced structural and stratigraphic interpretations within 
the Niger Delta basin to investigate the potentiality of hydrocarbon deposits (Ochoma, et al, 2023). Ochoma, et al, 
(2023) used the variance edge analysis to delineate the prominent and subtle faults in the area. The sweetness values 


range from 0 to 22,500. The high sweetness regions in the seismic data indicate high amplitude which indicates the 


ISSN; 2582-5267 OPEN@ ACCESS 


AJ B S R Asian Journal of Basic Science & Research 
JOURNAL Volume 5, Issue 2, Pages 83-96, April-June 2023 


presence of hydrocarbon-bearing sand units. The RMS amplitude values range from 0 to 13,000 in the reservoirs. 
The root mean square amplitude analysis also indicates the presence of hydrocarbon in seismic data. The relative 


acoustic impedance analysis was used for delineating lithology variation in the seismic sections. 


This study is taken from Fuba Field, Depobelt, Niger Delta, Nigeria. The ultimate deliverable of this study was 
application of 3-D seismic attributes analysis for hydrocarbon prospectivity of the area. The major components of 
this study are: (a) Well Correlation performed in order to determine the continuity of the reservoir sand across the 
field. (b) Seismic Interpretation which involves well-to-seismic tie, fault mapping, horizon mapping, time surface 
generation, depth conversion, generation of seismic attributes, seismic attributes analysis, prediction of reservoir 
properties from the seismic attributes and delineating hydrocarbon prospects using seismic attributes as direct 
hydrocarbon indicators (DHIs). This aids in giving more insight into application of 3-D seismic attributes for 


hydrocarbon prospectivity. 
“= 2.0. Location and Geology of the Study Area 


The proposed study area Fuba Field is located in the onshore Niger Delta region. Figure 1 shows the map of the 
Niger Delta oilfields showing the study area. The Niger Delta lies between latitudes 4° N and 6° N and longitudes 3° 
Eand 9° E (Whiteman, 1982). The Delta ranks as one of the major oil and gas provinces globally, with an estimated 
ultimate recovery of 40 billion barrels of oil and 40 trillion cubic feet of gas (Adegoke et al., 2017). The coastal 
sedimentary basin of Nigeria has been the scene of three depositional cycles (Short and Stauble, 1967). [14]. The 
first began with a marine incursion in the middle Cretaceous and was terminated by a mild folding phase in 
Santonian time. The second included the growth of a proto-Niger delta during the Late Cretaceous and ended in a 
major Paleocene marine transgression. The third cycle, from Eocene to Recent, marked the continuous growth of 


the main Niger delta. 


Figure 1. Map of Niger Delta Oilfields showing the location of Fuba Field 


A new threefold lithostratigraphic subdivision is introduced for the Niger delta subsurface, comprising an upper 
sandy Benin Formation, an intervening unit of alternating sandstone and shale named the Agbada Formation, and a 


lower shaly Akata Formation. These three units extend across the whole delta and each ranges in age from early 
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Tertiary to Recent. They are related to the present outcrops and environments of deposition. A separate member of 
the Benin Formation is recognized in the Port Harcourt area. It is Miocene-Recent in age with a minimum thickness 
of more than 6,000ft (1829m) and made up of continental sands and sandstones (>90%) with few shale 
intercalations (Horsfall et al., 2017). Subsurface structures are described as resulting from movement under the 
influence of gravity and their distribution is related to growth stages of the delta (Ochoma et al., 2020). Rollover 
anticlines in front of growth faults form the main objectives of oil exploration, the hydrocarbons being found in 
sandstone reservoirs of the Agbada Formation. The oil in geological structures in the basin may be trapped in dip 


closures or against a synthetic or antithetic fault. 
“= 3.0. Methodology 


Well correlation involves lithologic description, picking top and base of sand-bodies, fluid discrimination and then 
linking these properties from one well to another based on similarity in trends. Correlation of reservoir sands was 
achieved using the top and base of reservoir sands picked. The correlation process was possible based on similarity 
in the behaviour of the gamma ray log. In the Niger Delta, the predominant lithologies are sands and shales. In 
order to discriminate between these two lithologies in the subsurface, the gamma ray log is used. After defining the 


lithologies, the resistivity log was used for discriminating the type of fluid occurring within the pores in the rocks. 


There are seven basic steps involved in seismic interpretation relevant to this study and they include well-to-seismic 
ties, fault Mapping, horizon mapping, time surface generation, velocity modelling, depth conversion and attributes 
generation. The sonic log, which is the reciprocal of velocity, was calibrated using the checkshot data. The 
calibration process is necessary in order to improve the quality of the sonic log because the sonic log is prone to 
washouts and other wellbore related issues. The results of calibrating the sonic log with the checkshot gives the 


calibrated sonic log. 


The calibrated sonic log is used along with the density log to generate an acoustic impedance (AJ) log. The acoustic 
impedance log is calculated for each layer of rock. The next step involves generating the reflectivity coefficient 
(RC) log. The RC is calculated and generated using the AI log. The RC log generated is then convolved with a 
wavelet to generate a synthetic seismogram which is comparable with the seismic data. The statistical wavelet 
utilized for convolution is extracted from the seismic data. The synthetic seismogram was generated. The 


mathematical expressions that govern the entire well-to-seismic tie workflow are presented below: 


pE=22 PM (2) 
Pr Vo t+ PV 
Synthetic Seismogram = P2 V2 PAM wavelet (3) 
P2V2t PY; 


where AI = acoustic impedance, RC = reflection coefficient, p = density; v = velocity. 


Faults were identified as discontinuities or breaks in the seismic reflections. Faults were mapped on both inline and 


cross-line directions. Horizons are continuous lateral reflection events that are truncated by fault lines. The horizon 
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interpretation process was conducted along both inline and crossline direction. At the end of the horizon mapping, a 
seed grid is generated which serves as an input for time surface generation. Time surfaces were generated using the 
seed grids gotten from the horizon mapping process. The third order polynomial velocity model was generated and 


used to depth convert the time surfaces of the reservoirs of interest. 


Petrel software platform 2017 was used for the 3D seismic interpretation, attribute visualization and well logs data 
analysis. Among the seismic attributes that have been used in the visualization of the geology of the subsurface are 
variance, root mean square amplitude, maximum amplitude average magnitude and maximum magnitude. The 
variance attribute analysis was applied to the seismic inline 8515 while root mean square amplitude, maximum 


amplitude, average magnitude and maximum magnitude were generated as surface attributes. 
3.1. Determination of Root Mean Square (RMS) Amplitude 


Root mean square (RMS) amplitude is used to obtain a scaled estimate of seismic trace envelope. It is obtained in 
the software by sliding a tapered window of N samples as the square root of the sum of all the trace value x squared. 


The RMS attribute computation in Petrel software makes use of the inbuilt formula: 


1 2 
a a (4) 
where X,,,. =root mean square amplitude, w, = window values, N = number of samples in the window, x = 
trace value. 
3.2. Variance (Edge Detection) Method 


In the Petrel software, the variance attribute uses an algorithm that computes the local variance of the seismic data 
through a multi-trace window with user-defined size. The local variance is computed from horizontal sub-slices for 
each voxel. A vertical window was used for smoothing the computed variance and the observed amplitude 
normalized. The variance attribute measures the horizontal continuity of the amplitude that is the amplitude 


difference of the individual traces from their mean value within a gliding CMP window. 
2_1X 2 
o =—)> (x, +4,,) (5) 
n fi=l 


pene 2 : $ : 
Where o =standard deviation, 0” =variance, nm =the number of observations, S, = frequency, x, =the variable, 
X,, =mean of x; . 


3.3. Determination of Maximum Amplitude 


Amplitude is the deviation a wave from zero crossing. Maximum positive amplitude is referred to as peak value. 


The maximum amplitude attribute computation in Petrel software makes use of the inbuilt formula: 
y = Asin(at + @) (6) 


where A= Amplitude of the seismic wave, ¢ = Phase of the seismic wave, @ = Angular frequency of the seismic 


wave and f = time period. 
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3.4. Determination of Maximum Magnitude 


The maximum magnitude attribute computation in Petrel software makes use of the inbuilt formula: 
obs ' 
Mya, = Meee + [LF y ()] dn (7) 


where m,,,, =the maximum possible seismic magnitude, m?”° 


max 


=the largest observed magnitude of the seismic, 


m,,,, = the threshold of completeness, F, (7m) = the cumulative distribution function with which the magnitudes 


mu 
greater than or equal to the threshold of completeness are independent, identically distributed random values. n = 
all the seismic occurrences with a magnitude greater than or equal to the threshold of completeness, m=a 


measured magnitude of the seismic. 
3.5. Determination of Average Magnitude 


The average magnitude attribute computation in Petrel software makes use of the inbuilt formula: 


M, =DAu (8) 
2 
M = ge Mg It (9) 


M, =seismic moment, D =the average fault displacement, A =the total area of the fault surface, .s =average 
rigidity (with respect to shear forces) of the rocks in the fault, M,, =moment magnitude. 

“= 4.0. Results and Discussion 

4.1. Reservoir Identification, Correlation and Well-To-Seismic Ties 


The results for lithology and reservoir identification are presented in (Figure 2). A total of fifteen sand bodies (A, B, 
C, D, E, F, G, H, I, J, K, L, M, N and O) were identified and correlated across all seven wells in the field. Three 
reservoir sands were selected for the purpose of this study (M, N and O). The resistivity logs which reveals the 
presence of hydrocarbons were used to identify the hydrocarbon bearing sands. On Figure 2, the sands are coloured 
yellow while shales are grey in colour. Figure 3 shows the quality of Well-1 checkshot utilized for well-to-seismic 
tie. The results for well-to-seismic tie conducted on Fuba field using density log, sonic log and checkshot of Well-1 
is presented in Figure 4. An extended white 2 wavelet was used to give a near perfect match between the seismic 


and synthetic seismogram. 
4.2. Fault and Horizon Interpretation 


The results for the interpreted faults in Fuba field are presented in Figure 5 shows both synthetic and antithetic faults 
interpreted along seismic inlines. Faults are more visible along the inline direction because this direction reveals the 
true dip position of geologic structures. The variance time slice was used to validate the interpreted faults as seen on 
Figure 6. All interpreted faults are normal synthetic and antithetic faults. A total of thirty-six faults were interpreted 
across the entire seismic data. Of the 36 interpreted faults, only F1 (synthetic fault) and F4 (antithetic fault) faults 
are regional, running from the top to bottom across the field. Hence, these faults play significant roles in trap 


formation at the upper, middle and lower sections of the field. 
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The results for the interpreted seismic horizons (Horizon M, Horizon N and Horizon O) are also presented in Figure 
5. On these horizons, the fault polygons were generated and eliminated. The horizons were used as inputs for the 
generation of reservoir time surfaces. The reservoir time surfaces (M, N and O reservoirs) reveal that the reservoir 
structure is a collapsed crest, bounded by two regional faults (Fl and F4). Reservoir M, N and O time surfaces are 
truncated by two bounding faults and three minor inter-reservoir faults. The fault supported collapsed crestal 
structures are the closures identified on reservoir M, N and O reservoir surfaces respectively. The similarity in 
structure on reservoir M, N and O reveals that the field is structurally controlled by faults. The third order 
polynomial velocity function used for depth converting the studied surfaces is presented in Figure 7. The depth 
converted reservoir M, N and O surfaces are presented in Figures 8, 9 and 10 for the third order polynomial velocity 
function. The depth structure maps reveal that the reservoirs are anticlinal and fault supported. Reservoir M is found 
at a shallower depth from 10937 to 10997 ft, reservoir N is found at a depth ranging from 11213 to 11241 ft while 
reservoir O is found at a deeper depth ranging from 11681 to 11871 ft respectively. 
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Figure 2. Well Section Showing Reservoir Figure 3. Quality of Well-1 Checkshot Utilized 
Identified and Correlated Across Fuba Field for Well-to-seismic Tie 


Figure 4. Synthetic Seismogram Generation and 
Well-to-seismic Tie Conducted for Fuba Field 
using Well-1 Checkshot 


Figure 5. Faults and Horizons Interpreted 
Along Seismic Inline section 
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a Order Polynomial Function 
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Figure 7. Third Order Polynomial Velocity 
Model Utilized for Converting Reservoir 
Surfaces from Time to Depth 


Figure 6. Interpreted Faults Displayed on the 
Variance Time Slice 
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Figure 10. Reservoir Surface for (a) Time Surface 
O, and (b) Depth Surface O 


Figure 11. Variance Edge inline 8515 
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Figure 12. Generated Amplitude Maps for Reservoir Surface M (a) RMS Amplitude Map, (b) Maximum 
Amplitude Map, (c) Maximum Magnitude Map, and (d) Average Magnitude Map 
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Figure 13. Generated Amplitude Maps for Reservoir Surface N (a) RMS Amplitude Map, (b) Maximum 
Amplitude Map, (c) Maximum Magnitude Map, and (d) Average Magnitude Map 
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Figure 14. Generated Amplitude Maps for Reservoir Surface O (a) RMS Amplitude Map, (b) Maximum 
Amplitude Map, (c) Maximum Magnitude Map, and (d) Average Magnitude Map 
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Figure 15a. RMS Amplitude Map Showing Figure 15b. RMS Amplitude Map Showing Identified 
Identified Hydrocarbon Leads (A, B, C, D and E) Hydrocarbon Leads (A, B, C and D) away from the 


away from the Producing Zone “X” (for reservoir M) Producing Zone “X” (for reservoir N) 


A B CD EX 


Figure 15c. RMS Amplitude Map Showing Identified Hydrocarbon Leads (A, B, C, D and E) away from the 
Producing Zone “X” (for reservoir O) 


4.3. Seismic Attributes 


A series of seismic attributes such as variance edge, RMS amplitude, maximum amplitude, maximum magnitude 
and average magnitude were generated in Schlumberger’s Petrel® software interface to investigate potential 


structural and stratigraphic controls and also delineate hydrocarbon prospects within the study area. 


Figure 11 shows the computed variance attributes of the seismic section. The variance values range from 0.0 to 1.0. 
Values of variance equal to 1 represent discontinuities while a continuous seismic event is represented by the value 


of 0. The high values are denoted with red to yellow colorations. 


On the variance map, the areas dotted with blue, green, orange and pink colored lines signify values that correspond 
to the location of the discontinuity. The discontinuities may be interpreted as faults and boundaries as shown by the 
lines drawn on the variance attribute map (Law, and Chung, 2006). The variance edge enhanced the faults or 
sedimentological bodies within the seismic data volume. Furthermore, several bright spots are also delineated (in 
black circles and black ovals) which indicate high reflectivity sediments compared to their surroundings. These 
bright spots are an indication that a potential hydrocarbon trap might exist in the area. The darkest regions in the 
seismic section, which make vertical stripes, may be interpreted as faults or fractures. The zones with low variance 
values are due to similar seismic traces. Areas with red patches represent lineaments/discontinuities while grey 


areas represent the structural framework of the field. 
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The variance attribute is edge imaging and detection techniques. It is used for imaging discontinuity related to 
faulting or stratigraphy in seismic data. Variance attribute is proven to help in imaging of channels, fault zones, 


fractures, unconformities and the major sequence boundaries (Pigott, 2013). 


The RMS amplitude maps, maximum amplitude maps, maximum magnitude maps and average magnitude maps 
generated for all the studied surfaces (surface M, N and O) showed that strong amplitude anomalies exist at the 
central part of the field and parts of northeast, northwest, southwest and southeast sections (Figures 12, 13 and 14). 
The parts of the field where these strong amplitude anomalies were observed corresponds to the observed anticlinal 
structure in the field. Other parts of the field where the amplitude anomalies were observed also correspond to the 
identified structural closures found in the field. Results indicate that strong amplitude anomalies are structurally 
controlled as they are seen occurring near faults on structural maps. The RMS amplitude maps, maximum 
amplitude maps, maximum magnitude maps and average magnitude maps generated for all the studied surfaces 
presented were seen to have similar patterns of bright spot anomaly, which are likely associated with of facies 
and/or fluid content (Raef, et al, 2016). In addition, these anomalies are seen on structural high, which suggested 


prospect conformity with regional structural high. 


The results of the analysis of the RMS amplitude maps is presented in table 1. The RMS amplitude values range 
from 0 (purple) to 13,000 (red) in reservoir M, from 0 (purple) to 12,000 (red) in reservoir N and 0 (purple) to 9,000 
(red) in horizon O. The red-yellowish colour represents hydrocarbon sands. Some of these hydrocarbon sands were 
not detected in the original seismic section. The observed changes may be due to changes in lithology or fluid 
content. 


Table 1. RMS Amplitude Values for the Reservoir Surfaces 


Reservoir pune? enh 
Amplitude Values 
Reservoir-M 0 — 13,000 
Reservoir N 0 — 12,000 
Reservoir O 0 — 9,000 


The high values of RMS amplitudes may be related to high porous sands, which are potential hydrocarbon 
reservoirs. RMS amplitude is similar to reflection strength and it is used in seismic exploration for delineating 
bright spots and amplitude anomalies (Fozoa et al, 2018; Opara, and Osaki, 2018). The RMS amplitude is used for 
identifying coarser-grained facies, compaction related effects, and unconformities. The high values of RMS 
amplitudes circled (in orange circles) in the maps are interpreted as high porosity lithologies, such as porous sands. 
These high RMS amplitude segments are potential high quality hydrocarbon reservoirs. The high amplitude (in 
orange circles) in the seismic data conforms to the structures and confirms the presence of hydrocarbon (Ajisafe, 
and Ako, 2013). The high amplitude ranges from light blue to yellow and red coloration. Root mean square 


amplitude is used as a good indicator of the presence of hydrocarbon in seismic data. 
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4.4, Hydrocarbon Prospect Identification 


Evaluation of hydrocarbon prospectivity reveals that the study area has high fault density, which makes trapping 
structures available for hydrocarbon accumulation, thus providing possible hydrocarbon leads. Within the three 
reservoir intervals, the central part, northeast and southeast in an orange oval denoted with “X” has been previously 
evaluated and seen to have potential trap that turned into an oil/gas field when drilled, hence a good hydrocarbon 
prospect. This prospect is quite evidenced from the depth maps and several amplitude maps that were generated. 
Prospect “X” shows booming amplitude that is structurally controlled, which has given rise to several drilling 


campaigns within the zone (Figures 15(a)-15(c)). 


However, away from the “X” prospect this study has been able to identify other prospective zone characterized by 
good amplitude response and trapping structures. The root-mean-square amplitude attribute was used for the 
identification of prospective zones because it is more sensitive to direct hydrocarbon indicators (DHIs) than others 
and its pattern of bright spot anomaly superposition on the structural closures are discrete. Five hydrocarbon 
prospective zones each were identified for reservoirs M and O while four hydrocarbon prospective zones were 
identified for reservoir N (Figure 15 (a), (b) and (c). All the identified prospective zones in orange circles (A, B, C, 
D and E) in the reservoirs tops show two-way fault dependent closures. These been classified a hydrocarbon leads, 
since there are good trapping structures and hydrocarbon presence. Within reservoir M, zones A, B and C rank the 
best hydrocarbon lead considering the well-developed entrapment mechanisms for hydrocarbon accumulation and 
moderate to high amplitude response. Within reservoir N, zone C ranks the best hydrocarbon lead considering the 
well-developed entrapment mechanisms for hydrocarbon accumulation and moderate to high amplitude response. 
For reservoir O, zones B, D and E rank the best hydrocarbon lead considering the well-developed entrapment 
mechanisms with moderate to high amplitude response. This study has shown that, away from the producing zone 
(prospect “X,” which is common to all the three reservoir tops and characterized by anticlinal structure at the central 
part of the field), the application of seismic attribute analysis has unraveled several hydrocarbon prospective zones 


that could be further revalidated and evaluated to hydrocarbon prospects. 
“= 5.0. Conclusion 


A total of fifteen sand bodies (A, B, C, D, E, F, G, H, I, J, K, L, M N and O) were identified and correlated across all 
seven wells in the field. All interpreted faults are normal synthetic and antithetic faults. A total of thirty-six faults 
were interpreted across the entire seismic data. Of the thirty-six interpreted faults, only F1 (synthetic fault) and F4 
(antithetic fault) faults are regional, running from the top to bottom across the field. Hence, these faults play 
significant roles in trap formation at the upper, middle and lower sections of the field. Three horizons (M, N and O) 
were selected for the study. The seismic attributes interpreted include variance, root mean square amplitude, 
maximum amplitude, maximum magnitude and average magnitude. The variance values ranges from 0.0 to 1.0. 
The variance edge revealed the subtle structures and faults in the seismic section. The RMS amplitude values range 
from 9,000 to 13,000 in the reservoirs. The RMS amplitude, maximum amplitude, maximum magnitude and 
average magnitude results highlighted the hydrocarbon zones. The RMS, maximum amplitude, maximum 


magnitude and average magnitude attributes maps generated for reservoir tops of interest were seen to have similar 


ISSN: 2582-5267 OPEN@ ACCESS 


A J B S R Asian Journal of Basic Science & Research 
JOURNAL Volume 5, Issue 2, Pages 83-96, April-June 2023 


pattern of bright spot anomaly. These amplitude anomalies served as direct hydrocarbon indicators (DHIs), 
unravelling the presence and possible hydrocarbon prospective zones. Also, it indicates that away from the 
producing zone that there exist other hydrocarbon prospective zones. Several prospective hydrocarbon leads 
identified within three studied reservoirs (M, N, and O) intervals show good trapping structures such as anticlines 
which support hydrocarbon accumulation. Fault and horizon interpretation revealed that closures found on M, N 
and O reservoirs are collapsed crestal structures bounded by the two major faults. The depth structure maps reveal 
that the reservoirs are anticlinal and fault supported. Reservoir M is found at a shallower depth from 10937 to 10997 
ft, reservoir N is found at a depth ranging from 11213 to 11241 ft while reservoir O is found at a deeper depth 
ranging from 11681 to 11871 ft respectively. Extracted amplitude maps revealed that zones with good amplitude 
response, which indicates the presence hydrocarbon bearing reservoir section, are structurally controlled. The study 
has demonstrated that the application of seismic amplitude attribute, which offers insight on hydrocarbon presence 
and distribution of reservoir sand, could help unravel prospective intervals where there are no well controls. 
Hydrocarbon exploration and development risks can be reduced greatly with the outcome of seismic attributes 
extraction and analysis. It is recommended that careful quality control of the updated velocities using existing 
horizons should be performed so that the resulting velocity maps can be used for accurate depth conversion. Finally, 
fault seal analysis should be carried out to confirm that the suspected trapping faults are not leaking in which case 


they serve as conduits for hydrocarbon migrations rather than lateral barriers to hydrocarbon escape. 
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